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1. Introduction 

These lectures provide an introductory account of the theory of deep inelastic scat- 
tering. Many aspects of deep inelastic scattering can be analysed starting directly from 
QCD, without having to introduce phenomeno logical models. It is precisely these aspects 
of deep inelastic scattering that are most important, and are the focus of these lectures. 
These are a set of introductory lectures for students, and are not intended to be a review. 
I have therefore not given any references in the text. A few references are given at the end 
which might be useful as an introduction to the vast literature on this subject. 

The outline of these lectures is as follows. The kinematics of deep inelastic scattering is 
discussed in Sec. 2. Section 3 defines the leptonic tensor ', the hadronic tensor W^ u , and 
discusses the symmetry properties of and its decomposition into structure functions. 
The idea of scaling is covered in Sec. 4. The deep inelastic cross-section is derived in Sec. 5 
for spin-1/2 targets in terms of structure functions. A useful inequality on the structure 
function g\ that follows from unitarity is derived in Sec. 6. The physical interpretation of 
the structure functions in terms of virtual photon scattering is given in Sec. 7. This section 
also discusses the Compton amplitude T^ v , crossing symmetry, the optical theorem, and 
the analytic structure of the scattering amplitude. A brief discussion of deep inelastic 
scattering in the parton model can be found in Sec. 8. The operator product expansion 
and the derivation of the QCD sum rules for the moments of the structure functions are 
given in Sec. 9. This section is the most important part of these lectures, because it covers 
the properties of deep inelastic scattering that can be derived directly from QCD. Some 
applications of the results of Sec. 9, such as the Callan-Gross relation, the longitudinal 
structure function, the Ellis- Jaffe and Bjorken sum rules for gi, and the Wilczek-Wandzura 
contribution to gi are discussed in Sec. 10. Anomalous dimensions, and scaling violation 
are discussed in Sec. 11. Section 12 concludes with a brief summary of the structure 
functions b\{x) and A(x) for spin one targets. 

2. The Kinematics of Deep Inelastic Scattering 

In a typical deep inelastic scattering experiment, an incoming beam of leptons with 
energy E scatters off a fixed hadronic target. The energy and direction of the scattered 
lepton are measured in the detector, but the final hadronic state (usually denoted by X) 
is not measured experimentally. The lepton interacts with the hadron target through the 
exchange of a virtual photon; the target hadron absorbs the virtual photon, to produce the 



2 



final state X. If the target hadron remains intact, the process is elastic scattering. The 
deep inelastic region is where the target hadron is blown apart by the virtual photon, and 
fragments into many particles. I will only discuss the case of fixed-target deep inelastic 
scattering in detail. Deep inelastic scattering will soon be studied at HERA by colliding an 
electron beam with a proton beam. The kinematics for such colliding beam experiments 
is left as an exercise for the reader. 

The basic diagram for deep inelastic scattering is show schematically in fig. 1. There 
are numerous kinematic variables which are used in the discussion of deep inelastic scat- 
tering. In the definitions given below, I will pick the z axis to be along the incident lepton 
beam direction. (Warning: in later sections, I will pick the z axis to be along the direction 
of the virtual photon.) The kinematic variables are: 




P 

FIGURE 1. 

The basic diagram for deep inelastic lepton hadron scattering. The virtual photon mo- 
mentum is q. The final hadronic state is not measured, and is denoted by x. 

Kinematic Variables 

M The mass of the target hadron. The most important case is for a proton or neutron 
target, in which case M is the nucleon mass. 

E The energy of the incident lepton. 

k The momentum of the initial lepton. k = (E, 0, 0, E), if the lepton mass is neglected. 
O The solid angle into which the outgoing lepton is scattered. 
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E' The energy of the scattered lepton. 

k' The momentum of the scattered lepton, 

k! = (£", E' sin 6 cos 0, £" sin 6* sin 0, £" cos 6*) . 

p The momentum of the target, p = (M, 0, 0, 0), for a fixed target experiment. 

q = k — k', the momentum transfer in the scattering process, i.e. the momentum of the 
virtual photon. 

v = E — E' = p ■ q/M, the energy loss of the lepton. 

y = u/E = p ■ q/p ■ k, the fractional energy loss of the lepton. 

Q 2 = —q 2 = 2EE'(1 - cos 9) = AEE' sin 2 6/2. 

x = Q 2 /2Mv = Q 2 /2p ■ q = Q 2 /2MEy. 

u = 1/x. 

The variable x was first introduced by Bjorken, and is crucial to understanding deep 
inelastic scattering because we will see that QCD predicts that structure functions are 
functions of x and independent of Q 2 to leading order, a property known as scaling. Higher 
order corrections in QCD produce a small logarithmic Q 2 dependence of the structure 
functions which is calculable for large Q 2 since QCD is asymptotically free. I will therefore 
use the following definition of deep inelastic scattering: 

Deep inelastic scattering is the study of lepton-hadron scattering 
in the limit that x is Gxed, and Q 2 — > oo. 

The invariant mass of the final hadronic system X is 

M 2 X = (p + q) 2 = M 2 + 2p-q + q 2 . (2.1) 

The invariant mass of X must be at least that of a nucleon, since baryon number is 
conserved in the scattering process. This gives the inequality 

M 2 X >M 2 =>• M 2 + 2p ■ q - Q 2 > M 2 =>• x < 1. (2.2) 

Since Q 2 and v are both positive, x must also be positive. The lepton energy loss E — E' 
must be between zero and E, so the physically allowed kinematic region is 

0<x<l, 0<y<l. (2.3) 
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Eq. (2^2) can be written in the form 



x 



2p ■ q 



1 - 



M x — M 2 



2p ■ q 



(2.4) 



The value x = 1 implies that M x = M 2 , and so x = 1 corresponds to elastic scattering. 
Any fixed hadron state X with invariant mass M x contributes to the cross-section at the 
value of x obtained from Eq. (|2. J|). 



x 



1 + (M 2 - M 2 ) IQ 2 ' 



(2.5) 



In the deep inelastic limit Q 2 — > oo, so that any state X with fixed mass Mx gets driven 
to x = 1. In particular, all nucleon resonances such as the iV* get pushed to x = 1. 
The invariant mass of the hadronic state that contributes to the cross-section at a given 
value of Q 2 increases as x — > 0. For a fixed value of Q 2 , there is a small region around 
x = 1 of width Aq CD /Q 2 which probes hadronic resonances with masses around that of 
the nucleon. This resonance region is present in any real experiment since Q 2 is finite, 
but is not present in the formal deep inelastic limit. Outside the resonance region, the 
invariant mass M x is of order Q 2 . 

The experimental measurements give the cross-section as a function of the final lepton 
energy and scattering angle, d 2 a/dE'dQ. The results are often presented instead by giving 
the differential cross-section as a function of (x, Q 2 , 4>) or (x, y, <p). The Jacobian for 
converting between these cases is easily worked out using the definitions of the kinematic 
variables, 

1 

2MEy 2MEx 



d(x,Q 2 



d(x,y) 

d(x,y) 
d(£',cos#) 

Thus the cross-sections are related by 



O 2 

2MEx = — , 



y 



-2EE' 





d 2 



a 



dx dy i 



Mv\ d 2 o 



E' J dE'dQ 



E' 



0* 
y 



(2.6) 



d 2 a 



dx dQ^ 



(2.7) 



For most cases of interest, the cross-section will be independent of the azimuthal angle <f>. 
In these lectures, I will give the differential cross-section in (x, y, <fi). The variables x and 
y are Lorentz invariant, because they are defined using the inner products of four vectors. 
Thus the formula for the differential cross-section as a function of x, y, and <fi is also valid 
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for the kinematics at HERA, since cross-sections (and 0) are invariant under boosts along 
the collision axis. The only difference is that the transformation formula between (x, y), 
and {E' , 9) is modified. 

There are three independent variables which describe the kinematics, 9 and 
4>. The dependence on <fi is trivial, so I will only discuss the dependence on the other 
two independent variables. It is convenient to plot the allowed kinematic region in the 
(Q 2 /2ME) — {u/E)(= y) plane, as shown in fig. 2. The boundary of the physical region 
is given by the requirements that 

0<#<7T, 0<z/<£, 0<x<l. (2.8) 

Since x = Q 2 /2Mv = (Q 2 /2ME) / \v / 'E), the contours of constant x are straight lines 
through the origin with slope x. The relation between Q 2 and 9 can be written as 

Q 2 = 2EE'(l-cos9) =► _g_ = _L(E-^)(l-cos^). (2.9) 

Thus lines of constant 9 are straight lines passing through the point vjE = 1, and in- 
tersecting the Q 2 /2ME axis at Q 2 /2ME = (E/M)(l - cos9). Lines of fixed 9 become 
steeper as the beam energy increases, whereas lines of fixed x remain constant. 

The Q 2 dependence of the kinematic variables is crucial to understanding which terms 
are important in the deep inelastic limit. A generic point in the kinematic plane is given 
by some value of x and y. As E — > oo for a fixed value of x and y, the variables vjE 
and Q 2 /2ME are fixed. This implies that in the deep inelastic limit, v and E are each 
proportional to Q 2 , 

O 2 

v oc E oc (2.10) 
M v ' 

A generic point in the physical region has (1 — cos#) of order M/E, 

M M 2 

(l-cosfl) oc - oc ^, (2.11) 



so that the scattering angle 9 is of order M/Q, 

~Q 



M 

9 oc — , (2.12) 



and goes to zero as Q 2 — > oo. 

The allowed region in fig. 2 also shows that for fixed beam energy E, there is a limit 
to the Q 2 — x region which can be explored experimentally. The small x region is also the 



6 



G 2 



2ME 
A 



2E 
M 



E_ 
M 



\e=7i 
\ 

\ 

% 

\ 

\ 



\ n_7l 



x=7 



\ / 

y 




x= 



y 



v/E 







FIGURE 2. 

The kinematically allowed region for deep inelastic scattering is shown in grey. The dot- 
dashed lines are lines of constant scattering angle e, and the dashed lines are lines of 
constant x. In the deep inelastic limit, the intercept of the constant e lines with the 
vertical axis becomes infinite. 
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small Q 2 region, because lines of constant x approach the horizontal axis for small x. For 
a fixed value of x, the maximum allowed value of Q 2 is at the intersection of the line 9 = n 
with the line for fixed x. It is elementary to find the intersection point of the two lines, 



,2 / 2£ 



Qmax = 2M£x y 2E + Mx ) « 2M£x, (E » Mx). (2.13) 

To be in the deep inelastic region, one needs Q 2 to be larger than a few (GeV) 2 , so 
this places a limit on the smallest value of x accessible for a given beam energy. For 
example, with a 500 GeV lepton beam, and assuming Q 2 > 10 (GeV) 2 is large enough to 
be considered deep inelastic scattering, the smallest measurable value of x is 10 -2 . The 
factor of 2ME in Eq. ( |2.13 ) is proportional to the centre of mass energy squared, i.e. to 



s = (p + k) 2 . HERA will collide 28 GeV leptons on 820 GeV protons, so the value of 
s is 4(28) (820) (GeV) 2 . This will allow the measurement of structure functions down to 
x ~ 10 -4 . A similar measurement in a fixed target experiment would require a 50 TeV 
lepton beam. 

Finally, it is useful to have formulae for the different components of q as a function of 
x and y. The result is obtained by substituting for 9 in terms of Q 2 into the definition of 

q = k - k' = (E - E', E' sin 9 cos </>, E' sin 9 sin <fi,E — E' cos 9) , (2.14) 




3 Q 2 M 

q =2Mx-y +Mxy ' 



M 2 x 2 y 2 r 

-■Qy/l-y. 



(2.15) 



Q 2 (i-y) 

Note that q 1 - is of order Q, whereas q° and q 3 are both of order Q 2 /M. Also, q° and q 3 
are nearly equal, with q° — q 3 being of order M. These properties of Q will be useful when 
we compute cross-sections in the deep inelastic limit. 

3. Electroproduction, £^ u , and 

The basic Feynman graph for deep inelastic scattering is fig. 1. The scattering ampli- 
tude M. is given by 

iM = He) 2 ( Z |P) (tfltf (0) |*,^> (X\r h (0) \p, A) , (3.1) 



where sg is the polarisation of the initial lepton, A is the polarisation of the initial hadron, 
and jjf, j£ are the leptonic and hadronic electromagnetic currents respectively. The cou- 
pling constant e is not included in the definition of the current, and hence has been included 
explicitly in A4. The quantity A which describes the hadron polarisation can be chosen 
to be the value of the spin along an arbitrary quantisation axis (usually chosen to be the 
beam direction). For a spin 1/2 target, A = ±1/2. We will discuss the hadron polarisation 
in greater detail later in this section. 

The differential scattering cross-section is obtained from M. by squaring, and multi- 
plying by the phase space factors. The polarisation of the final lepton and hadron states 
are not measured, and have to be summed over. This gives the differential cross-section 

d 3 k' .. „,.„ IA1I 2 



E 



x 



{2E){2M){v reX = iy 
d 3 k' (2n) 4 5\k + p-k' - Px ) e 4 (3.2) 



[2n) 3 2E' (2E)(2M) Q 4 

x (p,A|#(0) \X) (X\f h (0) \p, A) (k,s e \j ilt (0) \k') (k'\j eu (0) \k,s £ ), 

where I have used the relation 

{a\fW = {(3\f\a), (3.3) 

which is true because the current is hermitian, y = j^. It is conventional to define the 
leptonic tensor by 

^= E (k'\jnO)\k,s e )(k,s e \j^O)\k'}, (3.4) 

final spin 



so that the leptonic current matrix elements in Eq. ( |3.2|) can be replaced by i^. The 
definition of the hadronic tensor is slightly more complicated. The tensor is defined 
by 

W^(p,q) x ,x = i- [ rf 4 xe^(p,A'|[^(x),r(0)]|p,A) , (3.5) 



4n _ 

where A and A' are the polarisations of the initial and final hadron. Inserting a complete 
set of states gives 

W^(p,q) x/X = ]-J2 [ d*xe i *- x \fa\ , \j»{x) \X) (X\f(0) \p, A) 

77 x J (3.6) 

-<p, A'|r(0) \X){X\j»{x) \p, A) 



9 



where the sum on X is a sum over the allowed phase space for the final state X. Translation 
invariance implies that 



(p, A'l^(x) \X) = (p, A'I^(O) \X)e^ p - px > x , 
(X\j»(x) |p, A) = (X\j»(p) \p, A) e^x-p)-*. 
Inserting Eq. (p.7|) into Eq. (|3.6|) gives 

W^(p, q)y X = ^ E Y 2 ^^ +P~Px) (p, A'| ^(0) \X) (X\f(0) \p, A) 

-(2n)W(q + Px -p) (p, X'\ f (0) \X) (X\j»(0) \p, A) 



(3.7) 



4-7T 



(3.f 



The only allowed final states are those with p° x > p°, since M\ > M 2 . Since q° > 0, only 
the first delta function in Eq. Q3.8f) can be satisfied, and the sum in reduces to the 
expression in Eq. (|3.2|) involving the hadronic currents, and the energy-momentum con- 
serving delta function (up to a factor of 1/4%). Only the first term in Eq. ( |3.6|) contributes, 
which is equivalent to the statement that one could have defined W^ v in Eq. ( |3.5|) simply 
as the matrix element of j' J '{x)j u (0) without the commutator. The reason for using the 
commutator is that then has a nicer analytic structure when continued away from 
the physical region which we will use extensively in Sec. 9. Substituting for and W^ v 
in Eq. (|3.2|) gives 

e 4 f d 3 k' 4tt l» v W^(p,k- k') xx 



da Q 4 J (2ti) 3 2E' (2E)(2M)(v Tel = l) ' (3 ' 9) 



so that 

f E' \ 



d 2 a 


e 4 


dE'dQ 


16tt 2 Q 4 


d 2 a 


e 4 


dx dy dtp 


16tt 2 Q 4 



(3.10) 



!/V W^(p,q) xx . 



All the information about the deep inelastic scattering cross-section is contained in the 
leptonic and hadronic tensors and W^ v . 

The Leptonic Tensor 

The leptonic tensor is trivial to compute since the leptons are pointlike fermions, 

t» v = u(k/)'y v u(k,8i)u{k,8i)'y> l u(k'). (3.11) 

final spin 
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The sum over final state spinors uses the identity 

u(k')u{k') = k" + m, (3.12) 

final spin 

where the spinors are normalised to 2E. 

The polarisation of a spin 1/2 particle can be described by a spin vector sg, defined 
in the rest frame of the particle by 

2st = u(k, s e ) au(k, s e ) = u(k, Si) 775 u(k, se), (3.13) 

where the states have the conventional relativistic norm of 2m in the rest frame. The 
vector s^ is defined in an arbitrary Lorentz frame by boosting si = (0, sg) from the rest 
frame, or equivalent ly, by 

2s£ = u(k, s*)7 M 75 u(k, s e ). (3.14) 

For a spin-1/2 particle at rest with spin up along the z axis, the spin vector is se = mz. 
This differs from the conventional normalisation of se by a factor of the fermion mass m. 
It is extremely useful to use this normalisation for the spin to avoid unnecessary factors 
of m appearing in the spin dependent cross-section. In the extreme relativistic limit, all 
mass effects can be neglected in the calculation, and all formulae should be written in 
terms of relativistic spinors normalised to 2E, without any additional factors of m. This 
suggests that one define spin by Eq. (|3.14| ) without any factors of m, and thus define a spin 



vector with the dimensions of mass, which is m times the conventional definition. With 
this normalisation for si, longitudinally polarised fermions in the extreme relativistic limit 
have S£ — Tiek, where k is the lepton momentum and Tie = ± is the lepton helicity. 

The initial state spinor product can be written in terms of the standard spin projection 
operator, 

u(k, s e ) u(k, s e ) = {jt + m) 1+ ^ m \ (3.15) 



Substituting Eq. fl3~T^ ) and flTI5D into Eq. ( |3~TT| ) gives 



= Tr(? + m e ) Y (it + rrn) 1 + 7 f /mi 

= 2 ( W + k u k'^ - gt* v (k ■ k' - mj) - i^ uaf3 q a Si P ) , ( 3 - 16 ) 
w 2 [k^k ,v + k u k'» - g^k ■ k' - ie^ a(i qci Sip) , 

neglecting the lepton mass. There are no factors of rrii in the spin term with our normali- 
sation convention Eq. (|3.14j ) . Note that the spin-independent part of the leptonic tensor is 
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symmetric in pi/, and the spin-dependent part is antisymmetric in pv. Thus an unpolarised 
lepton beam probes only the symmetric part of W^, and the polarisation asymmetry in 
the cross-sections probes only the antisymmetric part of W^ v . 

The Hadronic Tensor W^ u for Spin- 1/2 Targets 

The structure of the hadron relevant for deep inelastic scattering can be completely 
characterised by the hadronic tensor W^ v . Unlike the leptonic tensor which is known, W^ v 
cannot be computed directly from QCD because of non-perturbative effects in the strong 
interactions. The strong interactions are invariant under parity and time reversal, and 
these symmetries place restrictions on the form of W^. In this section, I will concentrate 
on hadronic targets with spin-1/2, which is the most important case experimentally. The 
tensor depends on the polarisations A, A', which for a spin-1/2 target can each take 
the values ±1/2. We will need to discuss the case where A / X' to describe an arbitrary 
target polarisation, even though it appears that only A = A' occurs in Eq. ( |3.1C1| ). 

The most general polarisation state of a spin J target in its rest frame can be written 
as a density matrix 

P= E E l A W(A'|. (3.17) 

A=-J A'=-J 

The states |A) and |A') transform as the spin-J representation under the rotation group. 
Thus the density matrix p transforms under the representation J®J = 0©1©...© 2 J. 
The density matrix of a spin-J target can thus be described by 2 J + 1 irreducible tensors. 
We will restrict ourselves in this section to J = 1/2, whose polarisation density matrix can 
be decomposed into irreducible tensors of spin and 1, i.e. by a scalar and a vector. It is 
easy to see that one can define a scalar and a vector from p by 

Pscalar = Tr p = 1. Rector = Tr p C = S h /M, (3.18) 

which defines Sh- p can be expressed in terms of these irreducible tensors, 

The only complication for targets with spin greater than 1/2 is that the decompositions 
analogous to Eqs. ( |3.18| ) and ( |3.19| ) are more complicated, which is why I have restricted 
the analysis of to spin-1/2 targets. A spin-1/2 polarised target can be described by 
an axial vector Sh, which has been defined with an additional factor of the target mass 
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M for the same reasons that we absorbed a factor of the leptonic mass into the lepton 
polarisation vector S£. The magnitude of Sh depends on the degree of polarisation of the 
target; for an unpolarised target, s*h = 0, for a 100% polarised target, \sh\ = M, and for 
a target with fractional polarisation /, \sh\ = fM. The four vector is defined to be 
(0, Sh) in the rest frame of the target, and is defined in other reference frames by a Lorentz 
boost. 

The tensor W^ u (p, q, Sh) for a spin- 1/2 target is defined by 

W^(p, q, s h ) = Px\> W^{p, q) yx = Trp (3.20) 

A, A' 

treating as a matrix in spin space with labels A, A'. Thus the most general tensor 
W^ v for a spin-1/2 target can be written as a function of the vectors p, q and Sh- Note 
that Sh contains the spin information of both the initial and final state spinors through 
the decomposition of the density matrix Eq. (|3.19| ). Thus Sh can only occur linearly in 
Wp,v\ there are no terms with two factors of Sh- 

The Decomposition o f W^ u into Structure Functions 

The deep inelastic structure functions are obtained by writing down the most general 
tensor decomposition of W^ v . The strong interactions are invariant under parity. The 
parity transform of the electromagnetic current is 

Vf(x)V~ l =r p (x P ), (3.21) 

with fjip = \x for fx = 0, and fxp = —fx for fx = 1,2, 3. To avoid writing complicated formulas, 
I will use the obvious convention j _1 = — j , etc., for the transformation of the various 
components of a vector under parity. The parity transform of the point x is xp = (t, —x), 
and the parity transform of the polarisation density matrix is 

VpV~ l = p P , Tr p P = 1, MTr p P a = s h , (3.22) 

since Sh is an axial vector. This implies 

W^(p,q,s h ) = ±-J d 4 xe^Trp[j»(x),f(0)} 

= ~L j d4xelQ ' XTrV P V ~ 1 WWifW?' 1 (3-23) 
= ^J d 4 xe^Trpp[j^(x P ),r-i 
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Defining qp = (q°, —q), (sh) p = (— s^), and using qp ■ xp = q ■ x, one can rewrite the 
last equality in Eq. (|3.23| ) as 

W» v {p,q,s h ) =W^(p P ,q P ,(s h )p), (3.24) 

which is the condition that is parity invariant. 

Invariance under time-reversal is only slightly more complicated, because time reversal 
is antiunitary. The time-reversal of momentum is defined by pp = (p°,—p), spin by 
( s h)T = ( s hi ~Sh)i an d coordinates by xp = (—t,x). Note that time reversal takes x° to 
— x , but still takes p° to +p°, i.e. time reversal changes the sign of time, but not of energy. 
Time reversal invariance then implies that 

W^( Pl q,s h ) = i- J d A xe^Trp[j»(x),f(0)} 

= ^J d A xe^ (Tr TpT- 1 T[f{x),f (O)]^ 1 )* (3.25) 

= ±-J d'xe^ (Trp T [j^(x T ),fT(0)])\ 

where the complex conjugation in the second line is present because time reversal is an- 
tiunitary. Using the identity q ■ x = —qp ■ xp, we get 



e 



x = (e iqT - XT )\ (3.26) 



so that time reversal invariance implies that 

W^(p, q, s h ) = W^ T (pp, q T , ( Sh ) T )*. (3.27) 
The electromagnetic currents and density matrix are hermitian, which implies 
W^(p,q,s h )* = ±- 1 rf 4 xe-^Trpt[^( x ) 5r( o)]t 

= ^-j d 4 xe-^Trp[f(0),r(x)} (3.28) 
= i-| d i xe- i ^Trp[f(-x),r(Q)}=W^(p,q, Sh ), 



where the first line follows from the identity (Tr AS)* = Tr(AB)^ and the last line follows 
from translation invariance. Note that all three identities, Eqs. ( |3.24j ), (|3.27| ), and ( |3.28| ) 
relate W^ v for g° to a transformed W^ v at the same value of q°. 
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There is a final identity that relates for q° > to W^ v for q° < which follows 
from the antisymmetry of the commutator and translation invariance: 

W»"(p,q, s h ) = ~ j d 4 xe^Trp[f(0)J»(x)} 

= -^J d 4 xe^Trp[f(-x),r(0)} (3.29) 

This identity is the photon crossing symmetry identity for W^ v , since exchanging pv and 
letting q — > —q is equivalent to exchanging the outgoing and incoming photons. Finally 
current conservation, d fJl j fJ '(x) = 0, implies that 

q^W^ip, q, s h ) = q u W^(p, q, s h ) = 0, (3.30) 

which greatly restricts the possible tensor structures allowed for W^ v . 

The most general hadronic tensor for polarised deep inelastic scattering from spin- 
1/2 targets can be written in terms of p,q,Sh and the invariant tensors g a p and e a p\ a 
with the constraints Eqs. (|3.24| )-( |3.30j ). The algebra is straightforward but tedious. It is 
conventional to write in the form 



w i i^iv \ . *2 f p-q Qn\ ( p-qq v 

W >» = Fl {- 9 ^ + ^) + p—q{ P ^ ) \ Pv " 

H e^Ao-g s +- rrr e^Acr? {p-qs -s-qp), 

p ■ q (p ■ qy 



(3.31) 



where the coefficients of the different tensor structures are called structure functions. The 
leptonic current is also conserved, so that q^l^v = q v (.^v = 0. Thus it is convenient to 
simplify the expression for W^ u by omitting all and q u terms before contracting with 

Wfj, v = -Fi g^ H p^p v H t^ v \ a q s h 

v ' q V (3.32) 

+ -^-e Mi ,Acr? [P- qs h - s h - qp ) , 

(Alternatively, one can drop the q^ and g„ terms in l^ u . However, it is incorrect to drop 
the g M and q v terms in both and W^l) There are also structure functions W±, W 2 , G\ 
and G 2 proportional to F±, F 2 , g\ and g 2 which are sometimes used in the literature, 

MW 1 =F U uW 2 = F 2 , 

(3 33) 

M 2 vG\ = #1, Mv 2 G 2 = g 2 , 
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but I will never refer to these, because they do not scale in the deep inelastic limit. 

For a spin-1/2 target, the symmetric part of W^ u is independent of the hadron spin, 
and the spin dependent part of is antisymmetric in \xv. We have already seen that 
the symmetric part of is independent of the lepton spin, and the spin dependent part 
of l^v is antisymmetric in [iv. Thus the combination W^l^v has no terms which have 
only the hadron spin, or only the lepton spin; all terms contain either both or none. Thus 
the structure functions F\ and F2 can be measured using an unpolarised beam and target, 
but to measure the structure functions gi and gi requires both a polarised beam and a 
polarised target. There is no advantage to doing an experiment with only a polarised beam 
or only a polarised target. This is not true for target spins greater than 1/2. For example, 
for a spin one target, there is target polarisation dependence in the symmetric part of W Ml/ , 
which can be measured using a polarised target and an unpolarised beam. 

4. Scaling 

I have referred to scaling at several points so far. We are now in a position to under- 
stand what scaling means. The tensor defined in Eq. ( |3.5| ) is dimensionless, as are 
the structure functions defined in Eq. ( |3.32|) . The structure functions are dimensionless 
functions of the Lorentz invariant variables p 2 = M 2 , p ■ q, and q 2 . It is conventional 
to write them as functions of x = Q 2 /2p ■ q and Q 2 = —q 2 , so they can be written as 
dimensionless functions of the dimensionless variables x and Q 2 /M 2 . In elastic scattering 
there is a strong dependence on Q 2 /M 2 , and the elastic form factors fall off like a power 
of Q 2 /M 2 . It was thought that the same behaviour would persist for the deep inelastic 
structure functions. The scale M is a typical hadronic scale, at which confinement (and 
other non-perturbative) effects become important. Bjorken was the first to point out that 
if the constituents of the hadron were essentially free pointlike objects at high energies, 
then the hadronic scale M should be irrelevant, and the structure functions then only 
depend on x, and must be independent of Q 2 . This is the famous prediction of scaling. 

We now know that QCD is an asymptotically free theory, and the strong interaction 
coupling constant a s becomes small at short distances. Thus at large Q, non-perturbative 
effects (such as the hadronic mass scale) are irrelevant, and QCD is described by a dimen- 
sionless coupling constant. Thus one recovers the prediction of scaling in QCD, since there 
is no dimensionful scale in the problem. As is well known, scale invariance in a quantum 
field theory is broken because of quantum corrections. This introduces scaling violations 
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because of anomalous dimensions and the running of a s . Since a s is small at high ener- 
gies, these scaling violations can be reliably computed in QCD perturbation theory. Thus 
QCD predicts scaling violations (i.e. Q 2 dependence) for the structure functions which are 
calculable. We will discuss this in more detail in Sec. 11. 



5. The Cross-Section for Spin-1/2 Targets 



The cross-section for spin-1/2 targets can now be obtained by combining the expres- 
sions Eq. ( |3.32| ) for W^, Eq. ( |3.16[ ) for £^ v , Eq. fl3.1U| ) for the cross-section, and using the 
identity 

(5.1) 



fiua/3 



-2 ( 9 x9. 



a 
9a 9 X 



d 2 a 



dx dy ( 



e 4 ME 
4n 2 Q 4 



xy 2 F 1 + (l-y)F 2 +y 2 g 1 [2x 



.-> 2 i Bh-Si 
+ 92 



p ■ q p ■ q p ■ q 
p- s e q- s h 



(5.2) 



p ■ q p ■ q p ■ q 

As mentioned earlier, the spin dependent terms involve both si and Sh- Let us look at 
the cross-section for the two cases which are most interesting for current experiments, a 
longitudinally polarised lepton beam incident on a target which is either longitudinally 
or transversely polarised. For a longitudinally polarised lepton beam, the polarisation is 
se = 7~Ltk, where Tii = ± is the lepton helicity. The lepton polarisation terms in Eq. 
can be written as 



q ■ st 
p-q 



Tit 



q ■ k 
p-q 



-Tie x, 



P ■ se 
p-q 



Tit 



p ■ k 
p-q 



TU 

y ' 



Sh ■ S£ 
p-q 



Tit 



Sh ■ k 
p-q ' 



(5.3) 
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Longitudinally Polarised Target 

A target polarised along the incident beam direction has Sh = MTi^z, where Tih = ± 
for a target polarised parallel or antiparallel to the beam. s^, is dotted into either q 
or k both of which have and 3 components which are almost equal. Thus Sh can be 
replaced by — Tih P i n the evaluation of the cross-section in the deep inelastic limit. (The 
additional minus sign occurs because of the relative minus sign between the space and time 
components in a dot product such as k ■ s.) The result is 

^dy = [X2/2jPl + (1 " ^ ~ HiHhy{2 ~ y)X91 + ° ( M2 /^ 2 )] ' ^ 

where the azimuthal angle <fi has been integrated over, since the cross-section is independent 
of (p. (The same expression can be used for HERA kinematics with the replacement 
2ME — > s.) Thus the polarisation asymmetry in the cross-section can be used to measure 
the structure function g\. The effect of the structure function g 2 is suppressed relative to 
the leading terms in Eq. (|5.4|) by an additional factor of M 2 /Q 2 , and has been omitted. 



Transversely Polarised Target 

The polarisation vector of a transversely polarised target can be chosen to point along 
the x axis. In this case 



d 2 a e 4 ME 



dx dy dcj) 2nQ 4 



xy 2 F 1 + (1 - y)F 2 + 2x 2 HiH h ^V^V cos</> {g x y + 2g 2 ) 



(5.5) 



where I have used Eq. ( |2.15|) for the magnitude of the transverse component of q. The 
structure functions gi and g 2 are equally important for a transversely polarised target, 
and so an experiment with a transversely polarised target can be used to determine g 2 , 
once gi has been measured using a longitudinally polarised target. Note that the spin- 
dependent piece of the cross-section for a transversely polarised target is smaller than 
the spin-averaged piece by a factor of M/Q, unlike in Eq. (|5.4[ ) where both terms are 
comparable in magnitude. Thus a small misalignment of the target relative to the beam 
direction can lead to a large contamination of the cross-section asymmetry due to the 
fractional longitudinal polarisation. It is therefore best to measure the cos^> dependence 
of the polarisation asymmetry to measure g 2 using a transversely polarised target. 
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6. Structure Function Inequalities 



There are some simple inequalities among structure functions that follow from uni- 
tarity. The differential cross-section Eq. (|5.4|) must be positive over the entire kinematic 
region, and for any sign of Hh and Tit. This requires that 



xy 2 F 1 + (1 - y)F 2 > xy (2 - y) \ 9l \ > 0, 



for all values of x and y. In particular, the inequality at y = 1 implies that 



Fi > \gi\>0. 



The n th moment of a function / is defined by 



M n (f) = I dx x n ~ l f{x) 

>0 



(6.1) 



(6.2) 



(6.3) 



We will derive sum rules for the moments of F% and gi in Sec. 9. There are useful inequal- 
ities amongst the moments that must be satisfied because of the restriction, Eq. ( |6.2[ ). For 
example, 



\M n ( gi 



dx x n gi(x) 



< dx x n \gi{x) 
'o 



< / dxx n - l F 1 (x)=M n (F 1 ] 
Jo 



(6.4) 



using the triangle inequality \a + 6| < |o| + \b\. Thus any moment of g\ is bounded by the 
corresponding moment of F\. There is also the inequality 



M n {Fx) = [ dx x^F^x) > f dx x n ~ 1 xF 1 (x) = M n+ i(Fi) 
Jo Jo 

since F\ > and < x < 1. Thus we find the sequence of inequalities 



(6.5) 



Mi (Fx) > M 2 (F 1 ) > Mz{F\) > M 4 (Fi) 



M X (F X )> \M x {g 



M 2 (F 1 ) >\M 2 ( gi )\, M 3 (Fi)> |M 3 (<7i 



(6.6) 



(6.7) 



The boxed quantities are the moments for which we will derive QCD sum rules in Sec. 9. 
The second moment of Fx, M 2 (Fi) bounds all the higher moments of F±, as well as all 
moments of g\ for which there are sum rules, except the first. This has an important 
experimental consequence which will be studied in Sec. 10, since the sum rule for the first 
moment of g\ is the Ellis-Jaffe sum rule. 
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7. Helicity Amplitudes: The Interpretation of Structure Functions 

The structure functions F±, F2, gi, and gi denned in Eq. ( |3.32| ) have simple physical 
interpretations. To understand them better, it is useful to look at the amplitude for forward 
Compton scattering off a hadron target (see fig. 3) 

(T^)y X = i [ d^xe^* ( Pl \'\T {j^x)j v {x))\ Pl \) ■ (7.1) 




FIGURE 3. 
The forward Compton scattering amplitude. 



The tensor T Ml/ (p, q,Sh) can be written in terms (T^x'x by analogy with Eq. ( |3.20| ), 

T^ v (p,q,s h ) = ^2p\x'T fl „(p : q) X 'x = Tr pT^. (7.2) 



A, A' 



The tensor Tp U (p,q,Sh) has the same symmetry properties as W^ u (p, q, Sh), Eqs. fl3.24p , 
([3.271) , (|3.28| ), and the current conservation constraint Eq. ( |3.30|) . Thus T^ v (p,q,Sh) can 
be expanded in the same tensors as W^ u . The structure functions for T fJiV (p, g, s/J will be 
denoted by a tilde to distinguish them for the corresponding structure functions for W^: 



- r , i ■ I^A , F 2 ( Pj_qqA ( P ■ q q„ 

H e^xaq s +- ^ e^uXaq [p-qs - s ■ qp ) . 

p ■ q (p ■ qY 



(7.3) 



The photon crossing symmetry relation for T, 



T^ u {p, q, s h ) = T vli (p, -q, s h ), 
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(7.4) 



differs by a sign from Eq. ( 3.29Q for W^ V1 because exchanging j M and j" inside a commutator 



produces a minus sign, whereas exchanging them inside a time-ordered product does not. 
Eq. (|7.4[) implies that F x is an even function of u>, and F 2 , gi, and g 2 are odd functions of 

F x (-uj) = F x (u), F 2 (-u) = -F 2 (u), 

g x (-u) = -gi(w), h(-u) = -02 (<«>)• 
The Compton amplitude for photon-antinucleon scattering can be obtained from the 
photon-nucleon scattering amplitude by fermion crossing symmetry, i.e. by replacing an 
outgoing hadron with momentum p and spin Sh by an incoming antihadron with momentum 
— p and spin Sh- Thus the Compton amplitude for antinucleon scattering is T^ v (—p, q, Sh). 
This implies 

ZTiantinucleon I , ,\ TT'nucleon/' , ,N ZTiantinucleon ( ,\ TT'nucleon/' , ,N 

~antinucleon ^^j^ ~nucleon^ ~antinucleon^^^ ~nucleon^ 

using Eq. ( f^lf) and noting that u — > — u if p — > —p. Eqs. (|7.5|) and ( |7.6|) imply that the 
antinucleon structure functions are equal to the corresponding nucleon structure functions, 

ZTantinucleon I , ,\ TT'nucleon/', ,\ Trantinucleon /, ,N TT'nucleon/', ,N 

~antinucleon ~nucleon^^^ ~antinucleon ~nucleon^^^ 

This relation could have been obtained trivially by applying charge conjugation to T Ml/ , 
since it converts the nucleon to an antinucleon. There are no additional signs because the 
product of two electromagnetic currents is even under charge conjugation. 

The optical theorem implies that twice the imaginary part of the forward scattering 
amplitude is the total cross-section. There are factors of i in the expansion of and 
Wfj, u , so what the optical theorem implies in our case is that twice the imaginary part of 
the T^ v structure functions is equal to corresponding W MJ/ structure functions times 4%. 
The extra factor of 4rr is due to the extra factor of l/4n in the definition of relative 
to T^ v . The analytic structure of the amplitudes F x , etc. is shown in fig. 4 as a function 
of a; = 1/x. 

There are cuts in the physical region 1 < \u\ < 00. The optical theorem implies that 

ImF x (u + ie) = 2irF x (u), ImF 2 (w + ie) = 2itF 2 (uj), 
Imr~/i(u; + ie) = 2ng x (uj), Img 2 (uj + it) = (<*>)■ 
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FIGURE 4. 

The analytic structure of t m „ in the complex u plane. The discontinuity across the cuts 
i<H<oo is related to vk m „. 

The discontinuity of T^ u across the right hand cut gives for nucleon targets. The 
antinucleon structure functions F±, F2, gi and gi are equal to the corresponding nucleon 
structure functions, because of Eq. ( J7.7| ), a result which also follows trivially from charge 
conjugation. We also see from Eq. ( |7.6| ) that the discontinuity across the left hand cut in 
the complex ui plane is the physical antinucleon structure function. 

The forward Compton amplitude can be expanded in terms of helicity amplitudes. 
Pick the z axis to be along the direction of the virtual photon, and choose this axis to 
quantise the angular momentum of the target and photon. The initial and photon spin 
components along z will be called h and h' , and the initial and final target spin components 
will be called H and H' . The helicity amplitude for forward Compton scattering (see fig. 5): 

7h + target H -> y h , + target^, (7.9) 

will be denoted by A(h, H; h' , H'). Angular momentum conservation along the z axis 
implies that 

h + H = h' + H'. (7.10) 
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h H H' h' 





y p p y 

FIGURE 5. 

The incoming photon scatters off the target proton which is at rest. The helicity amplitude 
for this process is denoted by A(h,H;h',H'). 

Reflecting in the plane of the scattering does not change any of the momenta, but reverses 
all the spins. Thus parity invariance implies that 

A(h, H- ti, H') = A(-h, -H; -ti, -H'). (7.11) 

Time reversal exchanges initial and final states, and reverses the directions of momenta 
and spin. Since the helicity of the photons and the target is defined as the component of 
spin along the direction of motion of the photon, helicity does not change sign under time 
reversal. Thus time reversal invariance implies that 

A(h, H; ti, H') = A(h', H'; h, H). (7.12) 

The virtual photon has momentum q = (q°, 0, 0, g 3 ), with three possible polarisations, 

e£ = --^(0,l,i,0), h = +l, 

e^ = -^=(0,l,-z,0), h=-l, (7.13) 
4 = ^(</ 3 , 0,0,(7°), h = 0, 

where e± are spacelike, and eo is timelike. The fourth polarisation proportional to does 
not contribute since the electromagnetic current is conserved. Denoting H = 1/2 by f, 
and H = —1/2 by { and using Eqs. ( |7.10j )-( |7.12| ), we see that for a spin-1/2 target there 
are only four independent helicity amplitudes which can be chosen to be 



-4(+,T;+,T), A(+,U+A), ^(0,T;0,T), ^(+,|;0,t). (7.14) 
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There are four independent structure functions for spin-1/2 targets because there are four 
independent helicity amplitudes. The helicity amplitudes can be computed in terms of T p!/ 
using 

A(h,H;h>,H , )=e£Ur i T^(s), (7.15) 

with 

s = Mu(H')au(H). (7.16) 

For spin-1/2 targets, -u(f) = anc ^ u {i) = ^l^' com P u ^ a ^i° n °f Eq. (|7.15| ) is 

simple for spin-1/2 targets in the deep inelastic limit. Using Eq. (|7.3[) and neglecting terms 
of order M/Q, we find 

A{+,V,+, J \) = Fi-g 1 , A{+ i U+,i) = Fi+g 1 , 

A(0, T; 0, T) = F 2 /2xF 1 = F L , A{+, |; 0, |) = (2x/2xM/Q) {g x + g 2 ). 



For example, the coefficient of F\ in Eq. ( |7.3| ) is —g^v Thus the coefficient of F\ in 
Eq. ( |7.17[ ) is +1 for A(+, T; +» T) an d ^(+>|j+>l) because e + is spacelike, is —1 for 
^4(0, |; 0,|) because eo is timelike, and is zero for A(+, |; 0, |) because e + and eo are 
orthogonal. Taking imaginary parts of the left hand side gives identical formulae with 
F — > 2nF, etc. Thus we see that 

2nF 1 = X - Im [A(+, T; +, T) + A(+, I; +, |)] , (7.18) 
is the cross-section for scattering a transverse photon off an unpolarised target, 

2tv 9i = -i Im [A(+, T; +, T) - A(+, |; +, |)] (7.19) 

is the spin asymmetry in the scattering cross-section for a transverse photon, and 

2nF L = lm^(0, T; 0, T) = lm^l(0, 1; 0, |) (7.20) 

is the cross-section for scattering longitudinally polarised photons off an unpolarised target. 
Fl is called the longitudinal structure function. The physical interpretation of g 2 is obscure, 
the combination gi +g 2 is proportional to the "single helicity flip" amplitude A(+, j; 0, f)- 
Note that the single helicity flip amplitude is suppressed by a factor of M/Q relative to 
the diagonal helicity amplitudes. 

The helicity amplitude formalism is useful in discussing targets with spin greater than 
1/2. It provides a simple way of counting the number of independent structure functions 
and understanding their physical significance. 
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8. The Parton Model 



The parton model provides a simple physical picture which helps in interpreting the 
structure functions of deep inelastic scattering. The hadron target is considered to be made 
up of a number of free on-shell partons (quarks and gluons) in the deep inelastic limit. 
The scattering cross-section can then be computed in terms of the incoherent scattering 
of free quarks and gluons. The hadronic tensor can be computed in terms of parton 
distribution functions. For example, the distribution function it-|-(£) is the probability 
to find an up quark in the hadron with polarisation parallel to that of the hadron, and 
momentum £p. It is conventional to go to an infinite momentum frame when discussing 
the parton model, and to claim that partons are non-interacting in this infinite momentum 
frame. This is incorrect. Partons can be treated as free because of asymptotic freedom, not 
because of the choice of a particular Lorentz frame. At high energies, the QCD coupling 
constant becomes weak, and so one can treat parton interactions using perturbation theory. 
The main reason for going to the infinite momentum frame is to neglect target mass effects. 
A parton with momentum fraction £ has invariant mass £ 2 M 2 , which is not zero, and 
depends on £. The infinite momentum frame hides this embarrassment, which is why it 
is used to discuss the parton model. Instead of using the infinite momentum frame, I will 
work in the rest frame of the target, and consistently neglect all target mass effects. 

To lowest order in the strong interactions, only quarks scatter off the virtual pho- 
ton, since gluons do not have electric charge. The contribution to from a single 
polarised quark can be computed from the Feynman graph of fig. 6. The computation of 
the scattering amplitude can be obtained from 




FIGURE 6. 



The contribution to w M „ from a single quark with momentum fraction £. 
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i^v, Eq. ( p.!6| ), by the replacements k — > £p, k' — > p', s# — ► s and g — > — g. and multiplying 
by the square of the quark charge Q. Comparing Eq. ( |3.16| ) and Eq. ( |3.8| ), we see that 
has an additional factor of integration over the final particle phase space. Thus we 

find 

^-i« 2 /(^l^ VKP + '- p,) (8,) 
x 2 [^p' v + & u p ,fM - g^ip ■ p' + ie^ a(3 q a s p ] . 

The only subtlety is the factor of l/£. The proton states in are normalised to 2p°, 
whereas the parton states in ( |8.1|) are normalised to 2p°£ in the standard relativistic 
normalisation convention. This affects the particle flux in the computation of the cross- 
section. The extra factor of l/£ converts the parton flux to the correct normalisation for 
the proton flux. The integral in Eq. (|8.1|) is easily done using the identity 



J w?^ = l w^^'-'')'! w?2Fl s \ ( + 2fl)- (8 ' 2) 

Since the partons are massless, sp = h^pp, where h is the helicity of the parton, so that 

= l^p' v + - gTZp ■ p' + iht e^q a p p ] 5{i - x), 

& q ( 8 - 3 ) 
- [2^ VP" " ^"fp " Q + ihZ e^q a pp] 5(£ - x), 



2£p-q 



where p' has been replaced by £p + g, and the g M and q v terms have been dropped. To 
compare with W^, we need Eq. Q3.32|) neglecting target mass effects. Thus the target spin 
Sh in Eq. ( |3.32| ) can be replaced by Hp, where 7i is the helicity of the target, and p is the 



target momentum. Comparing with Eq. ( |3.32j ), we see that the contribution of quarks to 
the structure functions is 

F 1 = ^-8{t-x), F 2 = Q 2 ^5^-x), 9l = ^hHS^-x), g 2 = 0. (8.4) 

The antiquark contribution can be computed similarly, and is identical to the quark contri- 
bution. The total structure function is obtained by integrating with the quark distribution 
functions, q+(£), the probability to find a quark in the nucleon with momentum fraction 
£ and helicity equal to that of the nucleon (JiH = +), and g-(£), the probability to find 
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a quark in the nucleon with momentum fraction £, and helicity opposite to that of the 
nucleon {KH = — ), and the antiquark distribution functions q±(C) : 

Q2 

F l {x,Q 2 )=^^- (q+(x) + q-(x) + q+(x) + q_(x)) , 

i 

F 2 (x,Q 2 ) = J2Q 2 i x (q+(x) + q_(x) +q + (x) + q_(x)) , 

i (8.5) 

Q2 

gi(x,Q 2 ) = ^2^- (q+(x) - q-(x) + q + (x) - , 

i 

g 2 ( Xl Q 2 )=0 1 

where the sum is over all quark flavours which can be considered light compared with 
Q 2 . The most important feature of this result is that the structure functions depend only 
on x, and are independent of Q 2 . We also see that the structure functions satisfy the 
Callan-Gross relation, 

F 2 (x) = 2xF 1 (x) => F L (x) = 0. (8.6) 

The Callan-Gross relation will get corrections at order a s in QCD. The structure function 
g 2 vanishes identically; there is no simple interpretation for g 2 in the parton model. 

We see from Eq. ( |8.5| ) that F\ is the probability to find a quark in the hadron with 
momentum fraction x, and gi is the difference in probabilities to find a quark in the 
hadron with momentum fraction x with spin parallel and antiparallel to the hadron. Both 
quantities weight each quark flavour by the square of the electric charge. The combination 
Fi(x) — g±(x) depends only on q-(x). The reason is that F\(x) — gi(x) is proportional to 
the imaginary part of the helicity amplitude A(+, |; +, |). There are two possible diagrams 
which contribute to the Compton scattering amplitude, one where the quark absorbs the 
initial photon, and then emits the final photon, and the crossed diagram where the final 
photon is emitted before the initial photon is absorbed. The crossed diagram has no 
imaginary part since the energy of the initial and final photons is positive, and q 2 < 0. 
The intermediate state when a quark with J z = m absorbs a helicity + photon is a virtual 
quark with J z = h + 1. A virtual quark only has J z = ±1/2, so only J z = —1/2 quarks 
can absorb a helicity + photon. This implies that only q~(x) contributes to F\(x) —gx(x), 
and similarly, only q+(x) contributes to F±(x) + gi(x). 

Isosyin Symmetry and the Gottfried Sum Rule 
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It is conventional to use the symbol q±(x) to refer to the parton distributions in a 
proton target. The distributions in a neutron target are related to those in a proton target 
by isospin invariance. Making a 180° rotation in isospin space exchanges p <-> n and u <-> d. 
Thus one finds 

proton (X) = d±, neutron \% ) 5 ^i, proton (x) = d±, neutron \X 1 5 

proton(x) = neutron^)} *^±, proton(^) = ^db, neutron j (8-7) 

■S±, proton (x) ^i, neutron 5 proton ^i, neutron (x) 5 

Isospin invariance is broken by light quark masses, and by electromagnetism. The light 
quark mass effects are of order m q / Aqcd, and the electromagnetic corrections are of order 
ck/47t, so one expects isospin violation at the few percent level. There is no large isospin 
breaking effect proportional to m u /md', the relevant quantity in QCD is the ratio of a light 
quark mass to the hadronic scale Aqcd- There is no reason whatsoever for the "quark sea 
to be isospin symmetric," i.e. for 

Wproton(lE) = rfproton(x), (8.8) 

because the proton state is not isospin symmetric. One simple way to see this is to note 
that the production of uu pairs in the proton must be different from the production of dd 
pairs in the proton, since there are two "net" u quarks in the proton, and only one "net" 
d quarks, and any u quarks produced must be antisymmetric with respect to the existing 
quarks. 

There is a (false) sum rule for the first moment of F^ p — F^ 71 , the difference in i*\ 
between protons and neutrons. This is the Gottfried sum rule which can be derived as 
follows: 

f F^(x) - Fr(x) = J 1 [| Hx) +u(x)] + i [d(x) +d(x)} + § [s(x) + s(x)]] pioton 

- [§ [u(x) +u(x)] + i [d{x) +d(x)] + ± [s(x) + s(x) 



neutron 



= ij Hx) - d(x) + u(x) - d(x)] pioton 
= I (Warning : This is false). 

(8.9) 

The second line follows from Eqs. (|8.7|) , and the last line follows from using Eq. (|8.8| ) to 
write 1 1 

I J [u(x)-d(x)+u(x)-d(x)] pioton = l ^ [u(x)-d(x)] = l (8.10) 

but since Eq. (|8.8| ) is incorrect, so is the Gottfried sum rule. There is another reason why 
the Gottfried sum rule is incorrect; it is an attempt to derive a sum rule for the "wrong" 
moment of F\ as we will discuss later. I repeat, the Gottfried sum rule is false in QCD. 
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9. The Operator Product Expansion and Sum Rules 

The parton model provides an interpretation for the deep inelastic structure functions. 
One can instead derive sum rules for the structure functions directly from QCD using the 
operator product expansion. These sum rules are extremely important because they do not 
depend on any model of hadronic structure, and provide a direct test of QCD. They can 
be derived using only some very general results from quantum field theory. The starting 
point in the derivation is the time-ordered product of two currents 

t^ = ij d 4 xe^ x T(j^x)j u (0)). (9.1) 

The nucleon matrix element of t^ v gives the familiar Compton amplitude, 

(V) VA = (p,A , |Vb,A). (9.2) 

The analytic structure of T Mi , and its relation to was discussed in Sec. 7. 

The key idea that permits the computation of T^ u in QCD is the operator product 
expansion. Consider the product of two local operators, 

O a (x) O b (0). (9.3) 

In the limit that x — > 0, the operators are at practically the same point. In this limit, the 
operator product can be written as an expansion in local operators, 

lim O a (x) O b (0) = Y j c ahk (x) O fc (0). (9.4) 

k 

The coefficient functions depend on the separation x. The left hand side is completely 
equivalent to the right hand side as long as one does not probe the operator product on 
distance scales which are small compared with the separation x. Thus one can replace the 
product O a (x)Ob(0) in the computation of matrix elements by the expansion Eq. (|9j|) 
where the coefficients c a bk(x) are independent of the matrix elements, provided that the 
external states have momentum components which are small compared with the separation 
x. In QCD, the coupling constant is small at short distances because of asymptotic freedom. 
Thus the coefficient functions can be computed in perturbation theory, since all non- 
perturbative effects occur at scales which are much larger than x, and thus do not affect 
the computation of the coefficient functions. 
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The momentum space version of the operator product expansion is for the product 

f^OMOM. (9-5) 

In the limit that q — ► oo, the Fourier transform in Eq. ( |9.1| ) forces x — > 0, and again the 
operator product can be expanded in terms of local operators with coefficient functions 
that depend on q, 

hm / d 4 xe^ x O a (x) O b (0) = Vc*( 9 ) O k (0). (9.6) 

q-HX> J ^ 

k 

This expansion is valid for all matrix elements, provided q is much larger than the charac- 
teristic momentum in any of the external states. 

We will use the Fourier transform version of the operator product expansion, Eq. ( |9.6| ) . 
The product of two electromagnetic currents in Eq. ( |9.1| ) can be expanded in terms of a 
sum of local operators multiplied by coefficients which are functions of q. This expansion 
will be valid for matrix elements in the target, Eq. ( |9.2|) , provided that q is much larger 
than the typical hadronic mass scale Aqcd- The local operators in the operator product 
expansion for QCD are quark and gluon operators with arbitrary dimension d and spin n. 
An operator with spin n and dimension d can be written as 

0£r"», (9-7) 

where Od, n is symmetric and traceless in \x\ . . . (JL n . The matrix element of Od, n in the 
hadron target is proportional to 

M d ~ n - 2 S\p^ ...p""], (9.8) 

for a vector operator, and to 

M d-n-2 S [ s Ml p M2 . , . j (9.9) 

for an axial operator. S acts on a tensor to project out the completely symmetric traceless 
component. The power of M follows from dimensional analysis, since a hadron state 
with the conventional relativistic normalisation has dimension minus one. The coefficient 
functions in the operator product expansion are functions only of q. Thus the free indices 
on the operator O must be either n,u, or be contracted with q a . Every index on O 
contracted with q a produces a factor of p ■ q (or s ■ q) which is of order Q 2 /M in the deep 
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inelastic limit. An index \i or v is contracted with i^, and produces a factor of p ■ k or 
p • k' (or s • k or s • k'), both of which are also of order Q 2 /M in the deep inelastic limit. 
In addition, since t^ v has dimension two, the coefficient of O must have dimension Q 2 ~ d 
in the operator product expansion. Thus the contribution of any operator O to W^ u £^ u is 
of order 



c. 



Mi •■■Mn d.n 



Q Q 
Q " Q P 



.p^ r 



(P • 9)" n 2-d M d-n-2 (9.10) 

where the twist t is defined as 

twist = t = d — n = dimension — spin. (9-H) 

The most important operators in the operator product expansion are those with the lowest 
possible twist. Twist two operators contribute a finite amount to the structure functions 
in the deep inelastic limit, twist three contributions are suppressed by M/Q, etc. The 
fundamental fields in QCD are quark and gluon fields, so the operators in the operator 
product expansion can be written in terms of quark fields ip, the gluon field strength 
and the covariant derivative We can make a table listing the basic objects, with their 
dimension and twist: 

1> Guv 



d 
s 
t 



3/2 2 1 
1/2 1 1 
1 1 



Any gauge invariant operator must contain at least two quark fields, or two gluon field 
strength tensors. Thus the lowest possible twist is two. A twist two operator has either 
two i/Vs or two G^s and an arbitrary number of covariant derivatives. The indices of the 
covariant derivatives are not contracted, because an operator such as D 2 has twist two, 
whereas D a D 13 has twist zero. 

The first step in doing an operator product expansion is to determine all the linearly 
independent operators that can occur. We have just seen that the leading operators are 
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twist two quark and gluon operators. The Lorentz structure of the quark operators must 
be either ip^y^tfj or '^/ ; 7 ^i 75' ^ / , m the limit that light quark masses can be neglected, because 
the operator does not change chirality. The conventional basis for twist two quark 
operators is: 

°v,a' tn = 2{^) S^f a r i D^...D^^ a ) , (9.12) 

0% a ^ = 2 (l) S^h^D^.-.D^^} , (9.13) 

where the index a = (u, d, s) runs over the light quark flavours. The twist two gluon 
operators are more complicated, and I will give one such operator tower: 

M» = __m s[G^ a D^...D^G aa ^] . (9.14) 

In these lectures, I will only compute the operator product expansion to lowest order, so 
I do not need the gluon operators. 

The second step in doing an operator product expansion is to determine the coefficient 
functions of the operators. The best way to do this to to evaluate enough on-shell matrix 
elements to determine all the coefficients. Since we have argued that the coefficients can be 
computed reliably in QCD perturbation theory, we will evaluate the coefficients by taking 
matrix elements in on-shell quark and gluon states. I will only illustrate the computation 
of the coefficients to lowest non-trivial order (i.e. (a s ) ) in this lecture. The generic term 
in the operator product expansion can be written as 

jj~c q O q + c g O g , (9.15) 

where q and g refer to quark and gluon operators. Taking the matrix element of both sides 
in a free quark state, gives 

(q\jj\q)~c q {q\O q \q) + c g (q\O g \q). (9.16) 

The electromagnetic current is a quark operator. Thus the left hand side is of order (a s )°. 
The matrix element (q\ O q \q) is also of order (a s )°, whereas the matrix element (q\ O g \q) 
is of order (a^) 1 since there are at least two gluons in O g , each of which contributes a 
factor of the QCD coupling constant g to the matrix element. Thus one can determine c q 
to leading order by taking the quark matrix element of both sides of the operator product 
expansion, neglecting the gluon operators. I will work for simplicity in a theory with a 
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p+q 




q 



p-q 



p, s 



p, s 



q 



(a) 



(b) 



FIGURE 7. 

The lowest order diagrams contributing to the quark matrix element of the product of two 
electromagnetic currents. 

single quark flavour with charge one, and generalise the result to an arbitrary number of 
flavours with arbitrary charges at the end. 

The quark matrix element of the left hand side of the operator product expansion, 
Eq. ( |9.15p , is given by the Feynman graphs in fig. 7, 



M» v = iu(p, s) 7""(P. s) + iu(p, s) 7 "z- ^ 



Yu{p, s). 



(9.17) 



(p + q) 2 ' ' (p-q) 2 

Note that there is an overall factor of % because we are computing % times the time ordered 
product in Eq. ( |9.1| ) and there are no factors of — ie at the vertex. The crossed diagram 
(second term) can be obtained by the replacement \i <-> u, q — > — q from the direct diagram 
(first term), so I will concentrate on simplifying the first term. Expanding the denominator 
gives 



(P + q) 2 = 2p-q + q 2 =q 2 ( 1 + 



2p-q 

q 2 



q 2 (l-u) 



(9.18) 



since p 2 = for an on-shell massless quark. The numerator can be simplified using the 7 
matrix identity 

(9.19) 



7a75, 
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u(p, s) 7 ^ + 4h u u(p, s) = u(p, s) (p + qYY + (p + qf^ - g» v + i) 

(9.20) 

+ ie^ aA (p + g) a 7A75j«(p, s). 

For an on-shell massless quark, 

^>u(p, s) = 0, u(p,s)j\u(p,s) = 2p x , s) 7A75 s) = 2hp\, (9.21) 

where h is the quark helicity. Thus the p 1 and e^ aA p Q 7A75 terms both give zero. (A 
note of caution: There have been many mistakes made by authors trying to distinguish 
u(p, s) 7a 75 u(p, s) p a from u(p, s) 7^75 u(p, s) p\. They are equal because of the equations 
of motion.) 

Combining the various terms, and using 



[l-cO-^Xy, (9.22) 



n=0 



gives 



2 00 



•^direct -^^"^ ^ + (P + ^ - 9^P '9 + ih e» vaX q.Px] • (9.23) 



n=0 

To complete the operator product expansion, we need the free quark matrix element of 
the right hand side of the operator product. The matrix element of the quark operators 
Eq. ( |9.12p and ( |9.13| ) in a free quark state is 

(p, s\ O^-^ \p,s) =p f * 1 ...p fln , (9.24) 

(p, s\ OJJ 1 "* |p, s) = hp^ .. (9.25) 

The factors of z and 2 in Eqs. ( |9.12j ) and (|9.13| ) were chosen so that no such factors appear 
in the matrix elements. 

I will determine the coefficient functions for the spin dependent term in the operator 
product expansion, and leave the computation of the spin independent terms as a home- 
work problem. The spin dependent term on the left hand side of the operator product 
is 

2 00 

M^ v] = ~-J2 un ihe^ aX q a p x + (// <-> v, -q, u -> -u), (9.26) 
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since ui is odd in q. The crossed diagram causes half the terms to cancel, so that the matrix 
element is 



A - . A J— ■ /— 9 \ n 

" n=0,2,4 " n=0,2,4 V 9 / 



n=0,2,4 ^ v ^ 



(9.27) 



n=0,2,4 ^ y ' 

where the dummy index A has been replaced by in the last line. Replacing n — > n + 1, 
and permuting the dummy indices fix . . . \i n gives 

n=l,3,5 ^ 9 ' 

The coefficient functions in the operator product depend only on q, and the matrix elements 
depend only on p. We have separated the operator product into pieces which depend only 
on q and only on p. By comparing with Eq. ( |9.25| ) , we can write Eq. ( |9.28|) as 

^[/H = ^ 2 q ie» va ^q a (0^- M "), (9.29) 

n=l,3,5 ^ _ 5 - 1 

so that 

£ [H = ^ 2 * ze^ Q ^ g^O^ 1 -^. (9.30) 

n=l,3,5 ^ _ ^ ' 

This is the operator product expansion for the spin dependent part of t^ v . We have 
omitted quark flavour types and quark charges in our analysis. The correct answer for 
Eq. ( p.30| ) has O a replaced by Qa®A,ai where Q a is the charge of flavour a in units of 
e. Odd spin axial vector operators are charge conjugation even, and even spin axial vector 
operators are charge conjugation odd. Since electroproduction is a charge conjugation 
even process, only the odd spin axial vector operators appear in Eq. (|9.30|) . I will leave the 



computation of the symmetric part of from Eq. Q9.23 ) as an exercise. The contribution 



can be written in terms of only even spin vector operators. (Even spin vector operators 
are charge conjugation even, and odd spin vector operators are charge conjugation odd.) 
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The most general operator product expansion at twist two can be written in the form 

, ( _ n | QfJ,Qu \ ■ ■ -Qfi n v^ 9 ^(1) ^Mi-.-M^ 

t>iw— 2.^ \ q 2 J (-q 2 ) n 2^ ZL/ j,™^j,v 



71=2,4,- ^ -x , K -x J j 

oo 



n=2,4,. 
oo 



2^ Ze M^A Ml 9 (-g2) n 2^ °J> 



; / - V- , . . KS ; I 

(9.31) 



A 



n=l,3,. 

oo 



, \^ ^ ■ ■ ■ (4) n/1 l/ /( ( /,.>/»;;...,. 

+ 2^ (Z^ 2^ 2 °i,n U j 



n=2,4,- 



J, A 



The sum in the various terms is over twist two quark and gluon operators. What we have 
just shown is that Cj„ at leading order is zero for the gluon operators, and is Q 2 . for 
the quark axial operator of flavour type a, as can be seen by comparing Eq. ( 9.31 ) with 
Eq. ( |9.3(J| ). A similar calculation for the spin-independent pieces shows that at leading 
order Cn and Cn^ are both equal to Q 2 , for quark operators of flavour type a, and all 
gluon operator coefficients are zero. (The terms of type have only gluon contributions, 
and so vanish to this order. They contribute to the structure function A(x) mentioned in 
Sec. 12.) 

To compute T MI/ , we need the hadronic matrix element of the operator product ex- 
pansion, so we need to evaluate the hadronic matrix elements of the operators Eqs. ( |9.12| ) 
and ( |9.13| ). The matrix elements of these operators in the target are not known. We can 
parameterise the matrix elements in terms of a known tensor structure times an unknown 
normalisation coefficient. The matrix elements will be considered only in a spin-1/2 tar- 
get. For arbitrary spin targets, the tensor structure of the matrix elements is much more 
complicated. The matrix elements of Oy and O a can be written as 

(p, s\ Oy 1 -^ \p, s) = V n p^ . . .p"», (9.32) 

(p, s\ O'X 1 -^ \p, s) = A n S [s^ 1 . . .p^} . (9.33) 

Using these matrix elements essentially undoes our free quark computation, except for 
factors of V n and A n , which were equal to one for the free quark matrix elements. Once 
again, let me concentrate on the term antisymmetric in \xv. Substituting Eq. ( ]9.33| ) into 
Eq. (|9.31|) , one finds 

T [H = V 2 Ci 3) ie" uafXl q a ' " - A n S [s^ 1 . ..p^} . (9.34) 
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Let us write 



(9.35) 



where 



= - [s^p^ 2 ...p^" + p^s^ 2 . . + . ..p^p^ 2 . ..s^"} - s^ 1 . . (9.36) 
n — 1 1 1 

_gMlpM2 pUn _| pMl ^2 _|_ _| pfJ-lp^2 gUn 



n n n 

The tensor 7£ Ml "" M ™ has no completely symmetric part, and so has spin n — 1 rather than 
spin n. Thus the contribution of 7£ to deep inelastic scattering is twist three, rather than 
twist two, even though it came from the matrix element of a twist two operator. This is 
the Wilczek-Wandzura contribution to #2> as will be discussed in Sec. 10. Omitting the 71 
term in Eq. (|9.35| ) for the moment, and substituting into Eq. ( |9.34j) , we find 

oo 

n=l,3,5 Q ( 937 ) 

= ie^ aX q a s x 

p-q 



using the definition of g\ from Eq. (|7.3|) . This gives us the power series expansion 



,n-l 



n=1.3,5 

(9-38) 



n=l,3,5 

The operator product expansion has allowed us to compute gi as a power series in u about 
uj = in QCD. The n th term in the power series is due to an operator of twist two and 
spin n. The radius of convergence of the power series is \u\ = 1, since that is the location 
of the first singularity in the complex uj plane. This is precisely where the physical region 
begins. 

The structure functions can be computed near u = 0, but not in the physical region 
1 < M < oo. To understand this better, it is useful to use light-cone coordinates, 

q ± = -L(q°±q s ). (9.39) 
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Then the deep inelastic limit, x fixed, Q 2 — > oo, is the limit q + fixed, q — > oo, because 
Q 2 = 2q + q~ , and 

1 o+ 

a; = - = (9.40) 

Thus deep inelastic scattering is not a short distance process; it probes the structure of the 
hadron along the light-cone. However, if we look at the unphysical region and let uj — > 0, 
then q + — > oo (and we already had — > oo). The expansion of the scattering amplitude 
around u = is a short distance expansion, which is why one can compute it in QCD. 

We can relate g\ in the unphysical region to its value in the physical region using 
contour integration. The coefficient 2CA A n can be extracted by the contour integral 

over a circle in the complex u plane around the origin, as shown in fig. 8. The contour of 
integration can then be deformed to the contour C. This gives 




FIGURE 8. 

The contours used to derive the moment sum rules. 
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We now have to assume that the contour at infinity does not contribute to the integral, 
and can be neglected. In that case, only the discontinuities across the cuts contribute to 
the integral over contour C. Using the optical theorem, Eq. 



gi {u + ie) -gx{(ji>- it) = 4ni g x (u) , (9.43) 

where gx(u) is the experimentally measurable structure function, and using Eq. ( |7.5| ) 

= -01 (w), (9.44) 

we get 

/OO 7 

9i(») J^, (9-45) 

where the 1 term is from the right hand cut, and the (— l) n term is from the left hand cut. 
Substituting u = 1/x, and using the definition of a moment, Eq. (|67|), gives 



2M n (gi) = 2 f dx x n - 1 g 1 (x) = C™A n , n odd. (9.46) 
Jo 



There is no sum rule for n even, because Eq. (|9.45|) reduces to the identity = 0. The sum 



rules for F\ and F2 are obtained in direct analogy to the derivation of the g\ sum rule. 
The result is that 

2M n (F 1 )=C^V n , neven, 

/Fo\ ,x (9-47) 

2Mn \i^ ) =c * Vn ' neven - 

The moment sum rules relate a quantity defined at high energy such as F\ to a low energy 
quantity, the zero momentum transfer matrix element of a local operator. One can derive 
sum rules in QCD only for the even moments of F\, and the odd moments of g±. There are 
no sum rules for the other moments. An attempt to derive a sum rule for such a "wrong 
moment" leads to the identity = 0, which is correct, but not very useful. This is due to 
the symmetry properties of the structure functions under charge conjugation (or crossing), 
Eq. ([7.4|). The Gottfried sum rule discussed in Sec. 8 is an attempt to derive a sum rule 
for the first moment of F±, and is not true in QCD. It tries to relate a charge conjugation 
even structure function F\ to a charge conjugation odd quantity, the number of u and d 
quark in the proton. 
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The only questionable assumption in the derivation of the sum rules is that the contour 
at infinity does not contribute to the integral. It should be possible to tell experimentally 
whether this is true. If the contour at infinity is important, then the sum rules will diverge 
near x = 0, which can be checked by measuring the small x behaviour of the structure 
functions. The sum rules become more convergent for the higher moments, so any problems 
with convergence will only occur for the lowest few moments. There is no evidence for any 
convergence problems for the sum rules in Eqs. (|9.46|) and (|9.47| ). 



10. Applications of the Moment Sum Rules 

The previous section was a long and involved discussion of the derivation of the mo- 
ment sum rules from QCD. In this section, I would like to discuss some of the consequences 
of the sum rules. 

Ft, m QCD 

The sum rules for F\ and F2 imply that the moments of the longitudinal structure 
function Fl = F2 /2x — F\ are given by 

2M n (F L ) = [C^ - V n , n even. (10.1) 

We have already seen that at lowest order, Cn = Cn\ This implies that Fl = at 
leading order in QCD. This is the Callan-Gross relation that we obtained previously in 
Sec. 8 using the parton model. Fl is more interesting in QCD than in the parton model. 
At order a s , Cn^ 7^ Cn\ so Fl 7^ 0, but is suppressed relative to F-y by a s (Q). The same 
operator matrix elements V n occur in Fl and Fi- The coefficients Clp and can be 
computed in QCD perturbation theory, so one can predict Fl by using the measured value 
of -Fi to determine V n . This prediction agrees with experiment. 

The Ellis-Jaffe Sum Rule 

Another interesting application is the sum rule for the first moment of g±, the Ellis- 
Jaffe sum rule, 

2 / g 1 (x)dx = C} 3) Ai. (10.2) 
Jo 

At one loop the coefficient is, 

c (3) = 1 _oOQ)_ 

7T 
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The quantity A\ is denned by the matrix element 

2s M Ai = (p, s 
This is conventionally written as 

A 1 = ^Au+^Ad+^As, (10.5) 

where Ag is defined by 

2s" Aq = (p, s| Q7^7 5 Q |p, s) . (10.6) 
In terms of the parton distribution functions introduced in Sec. 8, 

&q= [q+(x)-q-(x) + q + (x)-q_(x)]. (10.7) 
Jo 

The notation is slightly confusing because q+(x) — q~(x) is sometimes referred to as Aq(x). 

The Ellis- Jaffe sum rule relates the first moment of the g± structure function to the 
matrix elements of axial currents. This sum rule is special because there is no twist two 
spin-one gauge invariant gluon operator, so the right hand side of the sum rule only has 
the quark axial current operator, even at higher order in QCD. For ep scattering, the linear 
combination Au — Ad in the proton is equal to the axial decay constant qa in neutron 
beta decay by isospin symmetry 

(proton| w7 M 7 5 w — d^^d (proton) = (proton| u^^d (neutron) = 2gAS >1 . (10.8) 

The other non-singlet combination Au + Ad — 2As is determined by using 577(3) symmetry 
and fitting to the F and D coefficients in hyperon semileptonic decay, 

Aw + Ad-2As = 3F - D, g A = F + D. (10.9) 

The Ellis-Jaffe sum rule can thus be written in the form 

/ gi (x)dx= -p— [9F-D + 6As] =0.126 ±0.010 ±0.015, (10.10) 
Jo 18 

where the experimentally measured value for the sum rule is determined by the EMC 
collaboration (Ref. 6). 
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4_ 1- 1_ 

-^7^75^ ± -c?7 M 7 5 rf ± -57^758 



\p,s) 



(10.4) 



The values of F and D can be obtained by fitting to the hyperon semileptonic data. 
The standard results used in the literature are now out of date, because the experimental 
results have changed significantly. I will use the values 



F = 0.47 ±0.04, D = 0.81 ± 0.03, (10.11) 

from the paper by Jaffe and Manohar in Ref. 8. The EMC measurement, the values of F 
and D, and the Ellis-Jaffe sum rule determine all three Ag's, 

Au = 0.74 ±0.10, Ad = -0.54 ±0.10, As = -0.20 ± 0.11, (10.12) 

and 

A£ = Au + Ad+ As = 0.01 ±0.29. (10.13) 

Much of the recent interest in the spin structure of the proton is because of the non-zero 
value for As, and the small value for AE. 

The structure function inequality Eq. (|6.2| ) implies that s(x)+s(x) > \As(x) + As(x)\. 
(This can be proved by noting that the inequality holds for arbitrary values of the electro- 
magnetic charges of the quarks, and in particular holds for Q u = Qd = 0, Q s ^ 0.) The 
corresponding moment inequalities, Eqs. (|6.6|) and fl6.7j ) imply that the first moment of 
s{x) + s(x), i.e. the momentum fraction carried by s quarks, bounds all other moments of 
s(x) + s(x) and As(x) +As(x), except the moment that appears in the Ellis-Jaffe sum rule. 
The strange quark momentum fraction has been measured experimentally to be 0.026, so 
all other strange quark matrix elements should be small. The value of As, Eq. ( |10.12|) , 
escapes this constraint. 

The Bjorken sum rule is obtained by taking the first moment of the difference gl P — gf n 
for proton and neutron targets. Using isospin invariance, we see that 



j\gt P (x)-gr(x))dx= 1 -C{ 3) 



4 1 . 7 1 

-Au + -Ad+ -As 
9 9 9 

J proton— neutron 

= lc[V(Au-A d)proton 
"6 Cl 9A - 

Any flavour singlet contribution such as As or gluons cancels out in the Bjorken sum rule. 
This sum rule will soon be tested experimentally. 
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The Wilczek-Wandzura contribution to q 2 

Let us return to the tensor 71 that we left out in the computation of g\ in Sec. 9. 
Inserting 1Z into the formula for T^, Eq. ( |9.34|) , 



n 

n=l,3,5 



(10.15) 

since each term with a ^2 • • • A*n index on s contributes to the second tensor structure. 
Comparing with the expansion Eq. (|7.3j ) of T^ v , we see that Eq. ( |10.15|) contributes to the 
structure function g 2 , 

00 



n=l,3,5 



1 

1 

n 



-01 M + / ) — 



(10.16) 



using Eq. ( |9.38| ) for g\. Taking the imaginary part of both sides converts gi t 2 into the 
experimentally measurable structure functions g± t 2- Since g± has an imaginary part only 
in the physical region, the integral over ui' can be restricted to the region 1 < ui' < uj. 
Changing the integration variable to x = 1/u gives the Wilczek-Wandzura relation 

f 1 dx' 

g 2 (x) = -gi(x) + gi(x')—. (10.17) 

J X x 

This contribution to g<i came from twist two operators, but we have already seen that it 
is really a twist three contribution because 1Z has spin n — 1 instead of n. There are also 
twist three operators in the operator product expansion that contribute to g 2l so the total 
contribution to g 2 can be written as 

92(x) = [#2(z)] Wilczek _ Wandzura + MaO]^ three ■ ( 10 - 18 ) 

Both terms are equally important, and there is no reason why the Wilczek-Wandzura piece 
should be taken as a prediction for g 2 . 

The Wilczek-Wandzura contribution arises because the twist two operator Eq. ( 9.13 ) 



is symmetric in all its indices. This weights the p ■ s term in T Ml/ by a factor of (n — l)/n 
relative to the q ■ s term, where n is the spin of the operator. The standard decomposition 
of T^ u , Eq. (|7.3| ) weights the p ■ s and q ■ s terms by coefficients which are independent 
of the power of uj in the structure function. It is this mismatch which leads to twist two 
operators producing an effectively twist three contribution to g 2 . 
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11. Anomalous Dimensions and Scaling Violation 

The Q 2 dependence of deep inelastic structure functions is an extremely important 
subject, but I do not have the time to analyse it in any detail in these lectures. I just want 
to give an outline of the QCD analysis. The starting point is the moment sum rules for 
the structure functions, 

M n (F(x))~c q O q + c g O g , (11.1) 

in terms of quark and gluon operators, where F stands for a generic structure function. 
In deep inelastic scattering, the renormalisation scale \i is chosen to be Q, to avoid large 
logarithms in the computation of the coefficient functions. Thus Eq. ( |1 1 . 1|) is more properly 
written as 

M n (F(x, Q 2 )) ~ c q (Q) O q (Q) + c g (Q) O g (Q). (11.2) 

The Q 2 evolution of the structure functions is due to the Q 2 dependence in the coefficients 
and operators in Eq. ( |11.2|) . The coefficient c q (Q 2 ) = l+0(a 8 (Q)) depends only weakly on 
Q 2 . The largest Q 2 dependence comes from operator renormalisation. The renormalisation 
group equations for the operator can be written in the form 

= __J_ 7wOMl ... Mn) (1L3) 

at one loop, where 7 n is a calculable anomalous dimension. Dividing by the coupling 
constant evolution equation 



d , 9 3 

— 9 = —b — 
a\x 167r 

and integrating, gives 



»— 9 = ~bit^, (11.4) 



0(Q) 



O(Q ), (11.5) 



_a s (Q )_ 

so we can relate the moment of the structure function at Qo to the moment at Q. This 
leads to a calculable Q 2 dependence of the structure functions. As Q 2 — > oo, 0(Q) — > 0, 
if 7n > 0. 

The actual situation is slightly more complicated. The quark operator can be divided 
into singlet and non-singlet pieces. The non-singlet pieces evolve as discussed above, but 
the singlet pieces can mix with gluon operators. In the singlet sector, one needs to solve 
a two dimensional matrix renormalisation group equation. The computation is similar to 
the one sketched above, and the moment evolution equation is a 2 x 2 matrix equation in 
terms of the eigenvectors and eigenvalues of the anomalous dimension matrix. 
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The Q 2 dependence of the Ellis-Jaffe sum rule is particularly simple. There is no 
gluon operator that can contribute to the sum rule, and thus there is no gluon operator 
which can mix with the singlet quark currents. The non-singlet currents are conserved, 
and have no anomalous dimension. The singlet current is not conserved because of the 
axial anomaly, and has an anomalous dimension at two loops, 

= (11 ' 6) 

where Nf is the number of light flavours. Integrating the two loop equation by dividing 
by Eq. (IIT1D , gives 



kinglet (Q) = ex P 



(a(Q) - a(Q )) 



islnglet(Qo). (H.7) 



As Q 2 — > oo, j^(Q) approaches a finite value in contrast to the case where the operator 
had a one loop anomalous dimension, Eq. (|11.5|) . The Q 2 variation due to the two loop 
anomalous dimension Eq. (|11.7| ) produces a very weak Q 2 dependence in E in Eq. ( |10.13|) . 

There are some constraints on the possible anomalous dimensions for the structure 
functions in QCD. The unitarity inequality Eq. (|6.2|) must be true at all Q 2 . This implies 
that the anomalous dimensions of the axial vector operators must be larger than the vector 
operators, so that for large Q 2 , the g± structure functions falls off at least as fast as F±. 
Left and right handed quarks have the same QCD couplings. This implies that non-singlet 
axial and vector operators have the same anomalous dimension, since they cannot mix 
with gluons. This is not true for the singlet operators, because the Dirac trace in a closed 
Fermi loop depends on whether or not there is an additional 75 at the vertex. However, at 
one loop, it is still possible to show that the qq element of the singlet anomalous dimension 
matrix for vector and axial operators is the same, since there are no closed Fermi loops in 
the Feynman graph for the anomalous dimension. 



12. Outlook and Conclusions 

There were several new experiments on spin dependent deep inelastic scattering dis- 
cussed at this Winter Institute. The g± structure function of the proton will be mea- 
sured more accurately, to check the EMC measurement. The g± structure function of the 
deuteron is also going to be measured. This, in combination with the EMC measurement, 
will give the first measurement of g± for the neutron, and should be able to test the Bjorken 
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sum rule. There are experiments at HERA to measure two new structure functions which 
exist for spin one targets, h\ (x ) and A(x), which I did not have time to discuss here. 
The new generation of experiments are now sensitive enough to test the structure of the 
proton that is not "trivial," i.e. effects from pair production such as the strange quark 
distribution, Fermi statistics effects such as u 7^ <i, etc. 

Let me briefly discuss the structure functions b\ and A for spin one targets. In terms 
of helicity amplitudes, the structure function b\ is proportion to 



and is the angular momentum two asymmetry in the cross-section for a transverse photon 
to scatter off a polarised spin one target. The structure function bi is a twist two structure 
function that can be measured with a polarised target and an unpolarised beam. (All even 
angular momentum structure functions can be measured this way; odd angular momentum 
structure functions such as gi ( J = 1) need both a polarised beam and polarised target.) 
The structure function A(x) for a spin one target is proportional to 



and changes the helicity of the incident photon and target by two units. It is also a twist 
two structure function, but gets no quark contribution because the helicity of a quark 
cannot be changed by two units. Thus the structure function is leading twist but order a s , 
because it arises from gluon operators in the terms of type cj 4 ^ in Eq. ( |9.31] ). A(x) also 
vanishes for higher spin targets that can be considered as a collection of non-interacting 
spin-1/2 nucleons, so it should provide some interesting information on the structure of 
higher spin targets. 

I have tried to avoid any mention of models in this talk, except for a brief discussion 
of the parton model in Sec. 9. Instead, I concentrated on those aspects of deep inelastic 
scattering that can be derived from QCD using only quantum field theory. If a model 
calculation disagrees with experiment it is the model which is wrong (or needs to be 
"improved"); but if a QCD field theory calculation disagrees with experiment, it is the 
experimental result which is wrong, or it is evidence for a new fundamental interaction. 

Finally, I would like to thank Professor F.C. Khanna for inviting me to the Lake 
Louise Winter Institute. This work was supported in part by Department of Energy grant 
DOE-FG03-90ER40546, and by a Presidential Young Investigator award from the National 
Science Foundation, PHY-8958081. 



&i oc [A(+,l;+l) + A(+,-l;+,-l)-2A(+,0;+0)], 



(12.1) 



A cx ^(+,-l;-,+l), 
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